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Context. During solar flares an enormous amount of energy is released, and the charged particles, like electrons, are accelerated. 
I/") , These non-thermal electrons interact with the plasma in various parts of solar flares, where the distribution function of electrons can 

therefore be non-Maxwellian. 

Aims. We focus on the non-thermal components of the electron distribution in the keV range and analyse high-energy resolution X-ray 
Ph I spectra detected by RESIK and RHESSI for three solar flares. 

■ Methods. In the 2-4 ke V range we assume that the electron distribution can be modelled by an n-distribution. Using a method of 

line-intensity ratios, we analyse allowed and satellite lines of Si observed by RESIK and estimate the parameters of this n-distribution. 
At higher energies we explore RHESSI bremsstrahlung spectra. Adopting a forward-fitting approach and thick-target approximation, 
we determine the characteristics of injected electron beams. 
Results. RHESSI non-thermal component associated with the electron beam is correlated well with presence of the non-thermal 
n-distribution obtained from the RESIK spectra. In addition, such an n-distribution occurs during radio bursts observed in the 0.61 - 
15.4 GHz range. Furthermore, we show that the n-distribution could also explain RHESSI emission below ~ 5 keV. Therefore, two 
independent diagnostics methods indicate the flare plasma being affected by the electron beam can have a non-thermal component in 
the ~ 2-5 keV range, which is described by the n-distribution well. Finally, spectral line analysis reveals that the n-distribution does 
not occupy the same location as the thermal component detected by RHESSI at ~ 10 keV. 

Key words. Sun: flares - Sun: X-rays, gamma rays - Sun: radio radiation - radiation mechanisms: non-thermal - 
methods: data analysis 

1. Introduction iTsallisI d!988l) and iLeubnerl d2002l) have shown that the k- 
~, , , ■ • , , distribution is a consequence of the generalised entropy favoured 
. The energy release during the solar flares is accompanied by fe non . ext ensive statistics, therefore having a deeper physical 
_ ; many phenomena, i.e. magnetic field ^connection and reconfig- meani Kasparova & Karlickv (2009) interpreted RHESSI X- 
, uration particle acceleration plasma heating, and electromag- a of flares usi ^-distributions. Th found 

- ■ netic radiation across the whole spectral range. that SQme j x sources c£m fee described b such a 
X-ray spectra in deka keV range allow us to diagnose a ther- distribution 
mal part of the flare plasma, as well as its non-thermal, high- 

energy tail. Recent high-energy resolution, hard X-ray spectra |Seely et alj (|1987|) analysed SOLFLEX line spectra of 
obtained by Reuve n Ramaty High -Energy Solar Spectroscopic Fexxv and Fexxiv during several flares. They found that en- 
Imager, RHESSI, dLin et al.ll2002h enable us to infer the cor- hanced intensities of the satellite lines can be explained by the 
responding electron spectrum. This can be achieved with a presence of non-thermal energy distributions that have higher 
model-independent approach, such as regularised techniques and narrower peaks than the Maxwell distribution. This type 
dPiana et alJl2003L iKontar et alJl2004h . Also, the forward-fitting of distribution functions has been used in the analysis of n on- 
method is often used after assuming that the spectrum of thermal distribution in laboratory plasma (|Hares et alj[l979|) as 
the electrons emitting in this energy range is the sum of an well. In this paper we refer to them as n-distributions. Recently, 
isothermal compo n ent an d a non-thermal power-la w distribu- |Dzifcakova & Karhckyl d2008|) have introduced a simple model 
tion: e.g., iBrownl Jl97lb and iHolman et al] d2003h . This ap- of a mono-energetic electron beam penetrating the solar atmo- 
proach is adopted here. sphere and creating a return current. They showed that adding 

In space plasma, the non-thermal electron distributions a drift velocity corresponding to the electric current results in a 
are commonly mod elled by ^-distributions, based on in sit u non-thermal electron distribution that can be approximated by 
measurements, e.g., IVandasi dl989l) ; iMaksimovic et all dl997l) . the n-distribution. Both the electron beam and the n-distribution 

influence the line intensities but each in a different way. The 

* Appendix A is available at electronic form at http://www.aanda.org latter enhances the intensities of satellite lines, while the elec- 
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Table 1. Spectral lines used for diagnostics of the non-thermal distribution from RESIK spectrometer. 



Line number 


Ion 


Wavelength [A] 


Transition 


Ratio 


1 


Si xiv 


5.22 


ls 2 Sl/2-3p 2 P 3 /2,l/2 


Si xiv / Si xin 


2 


Sixm 


5.68 


Is 2 l S -ls3p l Pj 


Si xin / Si xnd 


3 


Si xnd* 


5.82 


is 2 2p 2 P l/2 ,3/2 - U2 P (3P) 3p 2 3 /2,5/2 





Notes. ( * ) The satellite line excited by di-electronic recombination. 



tron beam affect s the in tensities of allowed lines dominantly. 
Dzifcakova et al. (2008) analysed spectra of a flare using this 
kind of non-thermal distribution. They showed that the spectra, 
calculated using the n-distribution with the parameters estimated 
from observations, could mimic the observed relative intensi- 
ties of allowed and satellite lines very well. Additionally, they 
showed that the high intensities of the satellite lines observed 
by RESIK cannot be reproduced by a multi-thermal distribution, 
which further supports the preferred choice of the n-distribution. 

In this paper we apply two independent diagnostic tech- 
niques to analyse soft and hard X-ray flare spectra. In the soft 
X-ray range (~2 keV), where the bulk of the electron distri- 
bution is important, we consider the n-distribution. In the hard 
X-ray range (above 6 keV), we assume a single power-law dis- 
tribution describing the high-energy tail. A Maxwell distribution 
accounts for a possible thermal component at keV energies. We 
analyse flare spectra detected by two spectrometers and deter- 
mine the parameters of these model distribution functions. 



t=10 7 K N e =10'° cm" 1 
1.2x10 7 [ ' ' ' ' 




0.1 1.0 10.0 

E [keV] 

Fig. 1. n-distribution functions for different values of the param- 
eter n=\ , 3, 5, 9, 13, 19. n — 1 represents the Maxwell distribu- 
tion. 



2. Model distribution functions 

In the present work we consider three types of distributions: a 
Maxwell, a single power law, and an n-distribution, and for their 
determination we use different diagnostic methods. The thermal 
component is generally present and dominates the X-ray emis- 
sion around 10 keV. We suppose it is an isothermal emission so 
we approximate the corresponding electron distribution function 
as Maxwell. 

We assume that the har d X-ray emiss ion is caused by the 
thick-target bremsstrahlung (lBrownlll97ll) of an electron beam 
injected into the plasma with a single power-law spectrum: 

F(E) = (S- D^(^) = C.E~ S , (1) 

where 6 is the spectral index, Eq the low-energy cutoff, Fj rep- 
resents the total rate of injected electrons of E > Eq, and C is 
a constant. Sy mbol F(E) in Eq. Q] is the electron rate spectrum 
(lBrownll2005l) in units [s _1 eV" 1 ]. 

At the lower energies of X-ray flare spectra, i.e. 2-2.5 
keV, we use a diagnostics based on soft X-ray line spectra 
without continuum. We follow t he diagnostic technique devel- 
oped by Dzifcakova et aQ d2008l) . which uses the normalized n- 
distribution (Fig.Q]) to describe the electron distribution: 

f n (E)dE = B a ^-EHkTy (:i - +l) exp(-E/kT)dE, (2) 

where B n = 7T^ 2 /(2F(^ + 1)) is the normalization constant, E 
the energy of free electrons, k Boltzmann constant, and n and 
T are the free parameters of the distribution. The parameter T 



in this distribution is not the thermodynamic temperature, al- 
though it is also given in Kelvin. The parameter n is a dimension- 
less number describes the degree of deviation from the Maxwell 
distribution function. If n=l, the distribution function becomes 
Maxwell. The mean energy of the n-distribution, Eq. [3] depends 
on two parameters n and T. To compare the non-thermal com- 
ponent with the thermal one, we use the pseudo-temperature, 
t, which relates the mean energy of the n- distribution to the 
mea n energy of the Maxwel l distribution (iDzifcakoval 119981: 
iDzifcakova & Kulinovall2001l) : 

E n = (n + 2)kT/2 = 3&T/2, (3) 
where t - (n + 2)773 . 

2.1. Diagnostics of the n-distribution 

The line ratios of the three spectral lines listed in Table [T] were 
used for the diagnostics of the n-distribution. The intensities of 
two allowed lines (Si xiv and Si xm) depend on excitation rates 
that are integrals of a collisional cross section, velocity, and 
an electron distribution function. The intensities of the allowed 
lines are thus sensitive to the high-energy tail part of the electron 
distribution function (Fig. [2j. 

Besides the excitation equilibrium, the intensity ratio of the 
allowed lines is affected by the ionization equilibrium as well. 
The electron ionization rate of the Si xm is also an integral quan- 
tity that depends on the shape of the electron distribution func- 
tion above the ionization energy of Si xm. However, the radia- 
tive recombination rate depends on the shape of the whole elec- 
tron distribution function. Consequently, the ratio of the allowed 
lines is mainly sensitive to the pseudo-temperature r but only 
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Fig. 2. Schematic figure demonstrates which parts of a distribu- 
tion function influence the intensities of spectral lines (Table [TJ 
used for diagnostics of the n-distribution. The hatched regions of 
the distribution influence the intensities of the allowed lines of 
ions Si xiv and Si xm. The discrete energy of the doubly excited 
state of Si xnd is marked by a thick vertical line. 
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Fig. 3. Synthetic line ratios are functions of both n and r. The 
plot shows the synthetic line ratios for parameter n = 1, 3, 5, 7, 
9, 1 1, and the curves of constant log r (grey numbers). The curve 
n — 1 represents the line ratios for the Maxwell distribution. 



weakly depends on the parameter n, see Fig. [3] The Si xiv / 
Si xm ratio is very sensitive to T or r, and it varies from ~ 0.1 
for log r=6.8 to ~ 2 for log r=7.2. The magnitude of changes 
in this range of r is smaller for the higher values of param- 
eter n. This behaviour corresponds to the changes in the ion- 
ization and excitat i on equilibrium because of the n-distribution 
(Dzifcakovl [19981: [Dzifcakova & T6thoval l2007l) . The satellite 
lines of Si xn (further as Si xnd) show different behaviours. Their 
intensities depend on atomic parameters and the number of elec- 
trons with the energy corresponding to the energy of the dou- 
bly excited state of ion Si xn above its ionization limit (Fig. [2]). 
Thus, the ratio Si xm/Si xnd is extremely sensitive to the shape of 
the electron distr i bution function. For detailed information see 
Dzifcakova et al] (1201 lah . who investigate the influence of the 
non-thermal plasma bulk and the high-energy tail on the Si ion- 
ization and its line spectra. The observed enhanced intensities of 



satellite lines (Dzifcakov a et alJ2008h require an increased num- 
ber of particles at Si xnd energy in comparison to the Maxwell 
distribution (Figs. Q] [2] and[3]). The n-distribution with higher n 
satisfies this condition. 

Although the electron distribution is non-thermal, this does 
not necessarily imply that the ionization is out of equilibrium. 
The ionization rates of Si XIV are of the order of 10~ 10 cm 3 s _1 , 
which results in the ionization time of about 1 s for the elec- 
tron density 10 10 cm~ 3 . As the results show (this paper and 
Dzifcakova et al. 2008), the typical time scale of changes in the 
derived plasma parameters is longer, several tens of seconds. 
Thus, the ionization equilibrium can be achieved if the condi- 
tions for formation of the non-thermal distributions change on 
a longer time scale than the ionization time. That seems to be 
the case for the present analysis, therefore the assumption of the 
ionization equilibrium is valid here. 

In our model the n-distribution represents the bulk of the 
plasma electron distribution. It does not include the enhanced 
high-energy tail, representing a beam of accelerated electrons, 
usually observed at hard X-ray energies. 

We would like to point that, although the /c-distribution also 
provides the enhanced number of particles in its high-energy tail, 
its shape in the region of 1 - 4 keV is not able to explain the 
enhanced intensities of the satellite lines. This can occur only 
when the gradient of the electron distribution function is high 
enough in this energy region, and this is not fulfilled either by 
/(•-distribution or the Maxwell one. 



3. Instruments and data analysis 

We analysed flare spectra observed by two X-ray spectrometers: 
REntgenovsky Spekt rometr s Izognutymi Kristalami - RESIK 
(ISvlwester et alJl2005l) and RHESSI. 



3.1. RESIK 

The first spectrometer, RESIK, is an uncollimated spectrome- 
ter with two bent quartz crystals that covers the 2.0-3.7 keV 
energy range. It has four spectral channels with nominal wave- 
length ranges: 3.4-3.8 A, 3.83-4.27 A, 4.35-4.86 A, and 
5.00-6.05 A. In the present work we analysed spectra within 
the 5. 15 -5.95 A spectral range coming from the fourth chan- 
nel where the allowed lines of ions S xv, Si xiv, Si xm, and the 
satellite lines of ion Sixn dominate. The raw RESIK spectra 
were reduced to the absolute flux units by incorp orating all the 
known instrumental factors dSvlwester et al.l2005l) . Then, the re- 
duced spectra were time-weighted according to the exposition 
times of individual spectra. To improve the signal-to-noise ra- 
tio, we averaged the spectra in nearly equidistant time inter- 
vals, see Tables IA.1I IA.3I and IA.5I The diagnostic method of 
RESIK data is based on ratios of line fluxes, see Sect. 12. 11 there- 
fore we subtracted a linear continuum from the time-weighted 
spectra. However, RESIK spectra are contaminated by an arti- 
ficial continuum flux that are like mounds. The exact origin of 
these 'mounds' is still uncle ar, but we removed them using the 
same method as described in Dzifcakova et al. (2008). The im- 
perfectly removed 'mounds' can lead to relatively high errors 
of Si xiv line intensity, in particular, thus causing higher uncer- 
tainty of n for the spectra registered during a decay phase when 
this line intensity is low. Finally, the line spectra were fitted by 
Gaussian profiles with a constant full width at half maximum 
(FWHM). The errors of the line fluxes from the count statistics 
and fitting of the spectral lines reach 35%. 



3 



Kulinova et al.: Diagnostics of non-thermal distributions in solar flares 



Diagnostics of the n-distribution from the RESIK spectra 
makes use of the preca lculated grid of synthetic Si line ratios 
(Dzifcakova et al. 2008). These synthetic line ratios are func- 
tions of both parameters n and t, see Fig. [3] Applying the spline 
interpolation to the synthetic line ratios, we determined the val- 
ues of n and log r with their errors from the measured line ratios. 

Our line ratio diagnostic method can reliably detect only 
n < 11. For higher values of n, the curves of constant n start to 
overlap and lead to ambiguity in n. Therefore for all the mea- 
sured line ratios lying below the curve n — 1 1 (Fig. [3), the pa- 
rameter n was set artificially to 1 1 and the upper limit of n to 
12. Such an upper limit was chosen only for display purposes, 
and physi cally it means th a t n > 1 1 . In contrast to the previous 
paper bv iDzifcakova et al l (120081) . log t was determined using 
both line ratios (Table[TJ. There is a lower limit for log t: Fig. [3] 
shows that only log t ^ 6.8 can be reliably detected. In the case 
of n = 11, log t and its upper limit were determined as the values 
corresponding to n = 11 and 12, respectively. 

Obviously, the n-distribution with any parameter n 1 1 means 
a deviation from the Maxwell distribution. However, the er- 
rors from the count statistics and from removing the 'mounds' 
mainly affect the intensities of weak lines, and the higher uncer- 
tainties can usually lead to overestimated values of n during the 
decay phase. Therefore, in this study we concentrate on the time 
intervals of n > 5. 

3.2. RHESSI 

The second spectrometer, RHESSI, was designed to investigate 
the flare spectra in the energy range from 3 keV up to 7 MeV. 
We were interested in non-thermal bremsstrahlung radiation pro- 
duced by an electron beam injected to ambient plasma. In all the 
analysed flares the enhanced hard X-ray emission was observed 
approximately up to 100 keV. Among the RHESSI flare spectra, 
we studied in detail those time intervals during which we could 
observe both thermal and non-thermal components. 

In the hard X-rays, the non-thermal components of the flares 
lasted from two to about four minutes. At later times the hard 
X-ray emission of the flares could be fitted by one or two ther- 
mal components. We assumed the emission is composed of con- 
tinuum emission caused by electron-proton bremsstrahlung plus 
contributi ons of Fe (~ 7 keV) and Fe/Ni (~ 8 keV) line com- 
plexes dPhillips II20041) . 

We fitted the spectra in the Object SPectral Executive 
(OSPEX) environment using the model for the isothermal emis- 
sion and lines below about 10 keV (vth function with the de- 
fault set of the solar coronal abundances, iFeldman et alJPl992h 
and a thick-target n on-thermal emissio n at higher energies (pho- 
ton Jhick function. iBrown et~a l. 2008) assuming a single power- 
law electron beam, Eq. d_Q. If necessary, the p He up jnod pseudo- 
functiorQ with fixed default parameters was used. For spectra 
without the attenuator, state AO, the drmjnod pseudo-functiorfl 
was applied. We analysed the output from a single RHESSI de- 
tector 4F only bec ause of its fine energy resolution of 0.98 keV 
dSmith et alj|2002l) . Also, it produces a photon spectrum close to 
the mean spectrum o f all front detectors usable for spectroscopy 
(Phillips et al. 2006). Moreover, the drmjnod function can only 
be used to fit single-detector spectra. The low-energy boundary 
for the fits was set according to the attenuator state: ~ 6 keV in 



1 http://prg.ssLberkeley.edu/~tohban/wild/index.php/Pileup_mod_- 
_Pseudo_function_for_correcting_pileup 

2 http://sprg.ssl.berkeley.edu/~tohban/wiki/index.php/Drm_mod_- 
_Pseudo_function_for_fine_tuning_RHESSI_DRM_parameters 



AO, ~ 8 keV in Al, and ~ 12 keV in A3. The upper boundary of 
the energy fitting ranges was set to the values where the flare data 
were nearly equal to the background. In this way we obtained the 
time evolution of temperature of the thermal component T and 
the parameters of the injected electron beam, 6, Eq, and Fj. 

3.3. Radio data 

As complimentary data we used radio observations provided by 
USAF Radio Solar Telescope Networ k (RSTN) and Ondfejov 
radiospectrographs (Jifi cka et al.l fl993>. RSTN network consists 
of four solar radio observatories that monitor radio emission 
from the Sun on eight discrete frequencies (in GHz): 0.245, 
0.410, 0.610, 1.415, 2.695, 4.995, 8.800, and 15.400. The so- 
lar radio flux at these frequencies is recorded each second and is 
given in SFU units. Two radiospectrographs RT4 (2.0-4.5 GHz) 
and RT5 (0.8-2.0 GHz) at Ondfejov, Czech Republic, observed 
the Sun with 0.1s time resolution and sample the frequency band 
in 256 frequency channels. 



4. Flare events 

RESIK detected several flares between August 2001 and April 
2003. From its catalogue of flares we picked up three flares 
that were simultaneously observed by RHESSI and RSTN. One 
of the flares was observed by radiospectrographs at Ondfejov 
Observatory. 

4. 1 . 2003 February 22 flare 

The first flare was observed on February 22, 2003 in the active 
region NOAA 10290. According to Solar Event Reports it was 
of GOES class C5.8 with the soft X-ray beginning at 09:24 UT, 
and with its maximum at 09:29 UT and end at 09:31 UT. 

The results from the RESIK line spectra and RHESSI analy- 
sis are shown in Fig. |4] The top left-hand panel shows the com- 
parison of the Si xiv / Si xra vs. Si xm / Si xnd synthetic line ra- 
tios with the measured ones. The significant deviation of elec- 
tron distribution from the thermal one occurred slightly before 
the GOES soft X-ray maximum and during the time of increas- 
ing RESIK X-ray flux (09:27:43-09:29:03 UT), see Fig.g](right 
column). Later on, n ~ 3, which suggests slow thermalization of 
the flare plasma. 

Figure [4] (bottom left) shows that the log r decreased mono- 
tonically from its highest value, log r~7.3 (1.5 keV) observed 
in the rise phase of the flare. The log T from RHESSI is higher 
by ~0.1 dex than log t. Tables IaTI and |A21 list the determined 
RESIK and RHESSI parameters, respectively. 

The RSTN radio flux at several frequencies is plotted in 
Fig __l The peaks of radio bursts are observed within one minute 
before the GOES maximum, i.e. at ~ 09:28 UT. The flare 
shows a single peak at 0.61 -2.69 GHz, whereas there are mul- 
tiple peaks at 4.99- 15.40 GHz frequencies. The observed radio 
bursts correlate well with n > 5 detected from RESIK and with 
low 6 derived from RHESSI. Furthermore, the bursts occurred 
simultaneously with a group of dm-type III bursts observed in 
the time interval 09:28:16-09:28:33 UT, as shown in the radio 
spectrum in Fig. [6] Some of these bursts have positive frequency 
drifts and other negative ones. This indicates the electron beams 
moving in both directions: downwards and upwards in the solar 
atmosphere. 
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Fig. 4. RESIK and RHESSI analysis results of the 2003 February 22 flare. Top left: the synthetic line ratios for n=l, 3, 5, 7, 9, 11 
(solid, dotted, dashed, dot-dashed, dot-dot-dot dashed, and long-dashed lines, respectively), and the measured line ratios (stars) with 
their error bars. The arrows follow the time evolution of the RESIK spectra (0 marks the first observed spectrum, see Table lA.iT i. Top 
right: the time evolution of the parameter n (with error bars, upward arrows indicate n > 12) of the n-distribution and the total flux 
(dashed line) from all four RESIK channels. Bottom left: the time evolution of log t (diamonds) from RESIK and log T (triangles) 
from the RHESSI thermal component. Bottom right: the time evolution of the spectral index 6 (triangles) in comparison with the 
parameter n (diamonds and dashed line). This panel zooms in on the time evolution of n depicted in the panel above. The points in 
all the plots are depicted for the middle time of each spectrum. 



4.2. 2003 January 7 flare 



This M4.9 flare occurred in the active region NOAA 10251, 
starting in GOES soft X-rays at 23:25 UT, reaching its maxi- 
mum at 23:33 UT and ending at 23:40 UT. During this flare, 
see Fig. [7] parameter n started to rise about five minutes before 
the GOES soft X-ray maximum and remained > 5 until the end 
of the soft X-ray flare, except for a single drop at ~ 23:29 UT. 
After 23:42 UT it decreased below five. The parameter n gener- 
ally follows the total RESIK flux plotted in the top right of the 
Fig. [7j This is different from the previous flare where the max- 
imum of the total RESIK flux occurred about two minutes after 
the maximum values of n. The spectral index S shows a distinct 
time evolution from two other analysed flares. Its value stays al- 
most constant, ~ 7, with two noticeable drops, i.e. local minima, 
at 23:30 UT and 23:32 UT. 



On the other hand, the radio data suggest correlations with 
the time behaviour of 6 and n. The radio burst observed by RSTN 
(Fig.0 displays two peaks at frequencies 2.69- 15.40 GHz ap- 
pearing before the soft X-ray maximum. The times of these ra- 
dio peaks, i.e. ~ 23:30 and ~ 23:32 UT, correspond well to the 
two drops in the spectral index 6. Furthermore, the parameter n 
shows an increase up to 11 at the time of the second RSTN peak 
and remains at that value for several minutes (see RESIK spec- 
tra 5 and 6 in Table lA~3l and Fig. [7j top right). Additionally, the 
enhanced radio flux is still visible up to 23:40 UT. During that 
time, RESIK spectra still indicate deviations from the Maxwell 
distribution, exhibiting n > 5 (Fig. [§). In contrast to it, RHESSI 
spectra suggest thermalization of the hard X-ray emitting plasma 
sooner, at ~ 23:33 UT, when 6 rises to ~ 9. Since then, no signif- 
icant non-thermal component is present in the RHESSI spectra. 

Time evolutions of log t and log T during the flare are plotted 
in the bottom left of the Fig. [7] The determined values of log t 
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Fig. 5. Radio bursts observed on February 22, 2003 at six fre- 
quencies 15.40, 8.80, 4.99, 2.69, 1.41, and 0.61 GHz by RSTN. 
The dashed vertical line at 09:29 UT marks the time of the GOES 
soft X-ray flare maximum. The heavy black line denotes the time 
interval with the parameter n > 5 determined from RESIK and 
the heavy grey line represents time intervals with n < 5. The 
short black line denotes the time interval when groups of dm- 
type III bursts were observed, see Fig. [6] 



span the range of 6.81-7.29 (0.6- 1.7 keV). After reaching its 
maximum at ~ 23:31 UT, i.e. two minutes before the soft X-ray 
maximum, t starts to decrease until the end of the flare. Again, 
log T obtained from RHESSI thermal component are higher than 
log t with values between 7.29-7.44 (1.7-2.4 keV). In contrast 
to r, T has a rather spiky character. Its maximum was observed 
before the soft X-ray maximum, too, at 23:32 UT. Tables lA3l 
and lA.4l list the determined RESIK and RHESSI parameters, re- 
spectively. 

4.3. 2002 October 4 flare 

The third flare was of GOES class M4.0 and occurred in the ac- 
tive region NOAA 10137. It started at 05:34 UT, reached the 
GOES maximum at 05:38 UT, and ended at 05:41 UT. For 
this flare the parameter n>5 was observed between ~ 05:35- 
05:42 UT (Fig. top right). After ~ 05:39 UT n started to de- 
crease. Like the 2003 January 7 flare, the global behaviour of n 
follows the RESIK total flux. On the other hand, the time evo- 
lution of the spectral index 5 (Fig. [9] bottom right) has a similar 
shape to the 2003 February 22 flare (Fig. |4] bottom right), but 




Fig. 6. Composite radio spectrum at 0.8-4.5 GHz observed on 
February 22, 2003 by the Ondfejov radiospectrographs show- 
ing a group of type III bursts in the time interval 09:28:16- 
09:28:33 UT. 



in this case the RHESSI spectra are much steeper, therefore 5 
has higher values. Fast variations of n and its large uncertainties 
during the rise and the maximum phase are not a real feature, 
but are probably due to low photon flux and related high statisti- 
cal errors. The time evolution of log T obtained from RHESSI 
is nearly constant, and its values are again higher than logr. 
Tables|X3]and|At]list the determined RESIK and RHESSI pa- 
rameters, respectively. 

Radio burst observed by RSTN (Fig.fTOt occurred before the 
flare's soft X-ray maximum and started at ~ 05:36 UT. The shape 
of the radio flux shows a double peak at frequencies 4.99 GHz 
and 8.80 GHz with the first peak occurring at about 05:36 UT 
and the second one about 05:37 UT. At the highest RSTN fre- 
quency, 15.40 GHz corresponding to the lower chromosphere, 
the burst shows a triple peak between 05:36-05:37 UT. In con- 
trast to this, the 2.69 GHz flux shows a single broad peak with 
a maximum at about 05:37 UT. There radio bursts correlate well 
with the presence of the RHESSI non-thermal component ob- 
served at 05:35:28-05:37:32 UT. Furthermore, enhanced radio 
flux at all four frequencies, as well as n-distribution of n > 5, 
were still visible after the GOES maximum. 

5. RHESSI and n-distribution 

If there is any non-thermal n-distribution really present in flar- 
ing plasma, its electrons should also contribute to the production 
of bremsstrahlung radiation. Therefore, we included the thin- 
target model of bremsstrahlung from n-distribution into fitting of 
RHESSI spectra and tried to test the parameters of n-distribution 
independently from RESIK. We performed such an analysis only 
for the first two flares, 2003 February 22 and 2003 January 7, as 
we needed to use the RHESSI spectra observed in AO attenua- 
tor state for this purpose. This allowed us to examine the lowest 
part of the RHESSI energy range that is closest to the energies 
observed by RESIK. Unlike analysis of the RESIK X-ray line 
spectra, we suppose that RHESSI will give us the information 
about continuum radiation produced by the n-distribution. 

We therefore revisited AO RHESSI fits, as seen in the first 
row in Tables IA.2I and IA.4I and Fig. IA.1I where the fits are 
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shown with their normalized residuals. Both fits started at 6 
keV and they had;^ 2 ~ 1 (x\ means the reduced^ 2 ). Then, us- 
ing the best-fit parameters as a starting point, we pushed the 
low-energy boundary down to 4 keV and let the parameters of 
the isothermal and the power-law component vary. (When either 
pileupjnod or drmjnod pseudo-functions were used, the values 
of their parameters were taken from the previous fits and kept 
fixed.) However, no satisfactory fits could be obtained in this 
way. In this way xll$ was found for both flares, and the residu- 
als were unacceptably large, showing a systematic pattern below 
lOkeV (Fig-E9. 

Therefore, in the last set of fits we added the thin-target 
bremsstrahlung radiation from the n-distribution. This com- 
ponent was calculated by a new fitting function, thinjidistr, 
which has the four parameters to fit: n, r, "emission mea- 
sure" (i.e. ambient proton density x volume x electron density 
in the n-distribution), and a high-energy cutoff. (The function, 
thinjidistr, was incorporated into SSW tree and can be accessed 
from OSPEX.) The high-energy cutoff for the n-distribution is 
introduced for a numerical purpose only. It was set to 1 MeV 
and fixed for all measurements. The initial values of other three 
parameters of the thinjidistr function were set to those obtained 
from the RESIK diagnostics. 



These fits (Fig. fTTI) improved significantly both in terms of 
residuals and;^ 2 . The component due to the n-distribution de- 
scribes the lowest energies of the fits, whereas the spectrum at 
higher energies is not affected, and its parameters are similar 
to those from the original fits starting from 6 keV. Tables IA.7I 
and IA. 81 summarize the parameters for individual flares and for 
all three fits: the original fit, the fit with the enlarged fitting range, 
and the one combined with the n-distribution bremsstrahlung. 
To complete the derived parameters, we also list RESIK results 
for the corresponding time intervals in those tables. Tables lA77l 
and lA.8l and Fig.[TT]show the fits for a specific case of the param- 
eter n being fixed at a value derived from RESIK. The reason is 
that the studied Si lines provide the best possibility for determin- 
ing the parameter n because the line ratios with the satellite lines 
are very sensitive to the shape of the distribution considered. In 
addition to the fits with fixed n, we also analysed fits setting n 
as a free parameter. Both sets of fits were similar in terms of 
X 1 and residuals. However, in the latter case the best fit was not 
well defined since a wide range of n-distribution parameters re- 
sulted in a similar value for xl- Such behaviour suggested that 
the parameter uncertainties could be large. 

In order to estimate the uncertainties of the fitted parame- 
ters, we examined the^ 2 values by keeping a selected parameter 
fixed within a range of values and by fitting the other param- 
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Fig. 8. 2003 January 7 flare. Notation is the same as in Fig. [5] 



eters (see also chUjnap method in OSPEX). Having done so, 
we found out that x 1 does not exhibit well localized minima for 
n-distribution parameters. Moreover, in some cases only their 
lower or upper limits could be estimated. This is a situation sim- 
ilar to the low-energy cutoff (and consequently the total flux) for 
which it is known that in some cases the spectral fits could only 
provide an upper limit estimate. A cutoff below the upper limit 
still fits the data because at low energies the non-thermal flu x 
is usually well below the thermal component dSui et al.ll2005l) . 
On the other hand, this x 1 mapping technique allowed us to ob- 
tain thermal parameters and spectral index 5 uncertainties, i.e. as 
values corresponding to / y^ lin + 1. However, we note that this ap- 
proach does not provide joint uncertainties of fitted parameters, 
but instead corresponds to a single parameter uncertainty inde- 
pendent of any knowledge about the others. To obtain a joint es- 
timate of several parameters at a 68% ("lcr") confidence level, a 
much larger region, y * . + Ay 2 , would have t o be investigated 
dLampton et all fl976; m 'N aka mura et al. Il2010h . But this is be- 
yond the scope of the paper. 

The estimations of the parameter uncertainties are given in 
Tables ]Kj\ and |A8] For the case of the 2003 February 22 
flare, the fit with n fixed at the RESIK value, i.e. n = 11, gives 
n-distribution parameters that are slightly above the estimated 
lcr range, i.e. n-2-9, Ict = (1.0- 1.5) keV (but well within 2<x). 
However, due to the weak dependence ofx 1 on the n-distribution 
parameters, such a fit could be considered acceptable. Although 
the uncertainties of n and t obtained either from RHESSI or 
RESIK data are large for both events, their ranges overlap. 



We also performed a test to check the reality of the excess 
emission detected in the lowest RHESSI energies. A spectrum 
from a different detector (3F) was analysed and resulting fitting 
parameters agreed well within the uncertainties with those dis- 
cussed above. Therefore, the excess emission, fitted here by the 
n-distribution, is not an instrumental artefact but a real feature. 

Using the reconstructed images, we estimated the electron 
densities in plasma for different electron distributions. From the 
RHESSI fits we have the information about the emission mea- 
sure for the thermal component (EM^) and for the n-distribution 
(EM n ), as well as the total rate of injected electrons (Ft) for the 
power law. The n-distribution appears at the lowest part of the 
RHESSI spectrum. We therefore reconstructed the images in the 
3-6 keV range for it, the 6-12 keV channel was used for the 
thermal component, and the 25 - 50 keV channel was dedicated 
to the hard X-ray emission. From the images we determined the 
area S within the 70% contour of the maximum and supposed 
that the particular volume of emitting plasma has a spherical 
shape with a circular cross section S . Both for the n-distribution 
and the Maxwell distribution we supposed that N p = N e and 
£Mth = (A^Vth and EM n = (N e n ) 2 V„, where index 'th' is for 
the thermal component and 'n' corresponds to the n-distribution. 
The 70% contours in the 3-6 keV and 6-12 keV channels were 
almost the same (Fig.[T2l. so we used only the 6-12 keV chan- 
nel for the electron density estimation, assuming V± = V n . The 
electron density for the power-law component was derived using 
the relation N^ L ~ Fj/(Shaid ^/2E c /m e ), where Shard is the area m 
the 25-50 keV hard X-ray channel. 

In both flares the electron densities for n-distribution are 
three to four times higher than the electron densities of the ther- 
mal component. Moreover, they are about three to four orders 
of magnitude higher than the electron densities that correspond 
to the power-law distribution. Furthermore, the thermal and n- 
distribution electron densities are similar in both events, while 
N e for the power-law distribution differs nearly by one order. 
The 2003 January 7 flare was a limb flare (Fig. [12] right), so we 
register only its projection. It is likely that the flare footpoints, 
from where the dominant contribution to the 25 - 50 keV chan- 
nel comes, are the most affected by that projection. 

6. Discussion 

Basically, we suppose that a non-thermality in soft X-ray range 
manifested by the n-distribution is connected with a non- 
thermality in the hard X-ray region. Then at a first sight, one 
could expect to see an 'X-pattern' in n and 6 plots, i.e. when 6 
is the lowest, n reaches the highest values. But, in reality this 
connection is not so simple and does not depend on relation be- 
tween n and 6 only. It very likely also depends on electron flux 
density (Dzifcakova & Karlickvi|20 08) if we assume that the n- 
distribution is produced because of the return current. None of 
the three analysed flares exhibits a clear 'X-pattern'. There is a 
hint of such tendency, i.e. a smooth decrease in n versus a rise 
in 6, for the simple compact flare from February 22, 2003 of the 
GOES class C (see Fig. [4] bottom right), however the large un- 
certainties of n do not allow us to consider this as significant. The 
other two M class flares were more complex. The 2002 October 
4 flare was a two ribbon flare occurring near the limb centre. The 
2003 January 7 flare was a limb flare exhibiting a rather com- 
plex magnetic configuration, possibl y with a plasmoid structure 
(Fig. [HI and iDzifcakova et alJl2008l Fig. 2). Additionally, both 
these flares show longer time occurrence of the non-thermal n- 
distribution than the power-law component observed in RHESSI 
spectra. 
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Table 2. Derived electron densities for Maxwell (A 7 *), power 
law (N^ L ), and n-distributions (N") 



Distribution 


February 22, 2003 


January 7, 2003 


Af [cm" 3 ] 


1.6xl0 10 


9.8xl0 9 


A£ [cm" 3 ] 


4.0xl0 10 


3.5X10 10 


jv c pl [cm- 3 ] 


4.7xl0 6 


2.8xl0 7 



In all three flares the time interval with n > 5 corresponds 
well with radio bursts observed at frequencies 4.99-8.80 GHz. 
Also, we found a good correlation of radio emission with hard 
X-rays related to electron beams. This behaviour is known and 
has be en st udied, e.g.. bv | Aschwanden et al.l (1 19951) . iBenz et a"D 
d2005l) . and lFarnfk & Karlickvl ( l2007l) . For the 2003 February 22 
flare, no enhanced radio emission was observed after the burst. 
On the other hand, the long duration enhanced radio emission 
was presented, together with n > 5 for both M class flares, but a 
correlation with 6 shows a shorter duration. 

For the 2003 February 22 flare, the type III bursts with posi- 
tive and negative frequency drifts were observed within the time 
interval with n>5. These bursts indicate the electron beams 
are moving downwards and upwards in the solar atmosphere. 
Because each electron beam needs to be a ccompanied also by 
the return current electrons (iKarli ckv 2009), such a radio spec- 
trum can be considered as evidence of the non-Maxwell electron 
distribution function representing an electron beam with a return 
current in the region of the radio source. Assuming radio emis- 
sion generated on the plasma frequency (800-4500 GHz, see 
Fig. |6j, these type III bursts thus show that deviations from the 
Maxwell distribution are present in the atmospheric layers with 
the plasma density N e from 7.9xl0 9 cm -3 to 2.5xlO n cm -3 , 
i.e. in the chromosphere through the transition layer to the low 
corona. The electron density of the n-distribution determined 
from the RHESSI spectral fits and images, see Table [2] agrees 
well with the above estimate from the radio emission. This fur- 
ther supports the idea of electron beams and of an accompany- 
ing return current that causes the n-distribution. Additionally, as- 
suming that the ratio A^ u /7V e = lO" 6 - lO" 7 ( Mel rosel Il980l 
p. 203 - 204), where A/' 11 is the plasma density of the electron 
beam generating type III bursts, then we can see that A 7 ' 11 is much 
lower, * 7.9 x 10 3 - 2.5 x 10 5 cm , than the plasma densities 
of the n-distributions presented in Table [2] Thus a contribution 
of the electron beam electrons producing type III bursts to the 
n-distribution is very small. 

Furthermore, we made the following comparison of pa- 
rameters of flares in February 22, 2003 and January 7, 2003. 
In the model in which the n-distribution is formed by elec- 
trons of the return current induced by the electron beam 
dDzifcakova & Karlickvl l2008h . we can write A^vj^. = Ar b v b , 
where A 7 " is the electron density corresponding to the n- 
distribution, Nb the density of the power-law tail (electron 
beam), v^ c the mean return current velocity in the proton co- 
ordinate system, and Vb the mean velocity of the power-law tail 
(beam). The n-distribution densities in both the flares are nearly 
the same (Table [2]). Therefore, if we assume the same velocity 
v^ c , i.e. the same n-distributions in both flares, then the mean 
velocity of the power-law tail (electron beam) vj in the 2003 
February 22 flare should be greater than that in the 2003 January 
7 flare. However, estimation of the mean velocities significantly 
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Fig. 10. 2002 October 4 flare. Notation is the same as in Fig. [5] 



depends on the low-energy cutoff that was determined with large 
uncertainty. Nevertheless, the power-law index 6 ~ 3 in the 2003 
February 22 flare is much smaller than that (6 ~ 6.5) in the 2003 
January 7 flare, which means that the mean velocity in the 2003 
February 22 flare is greater than in 2003 January 7, as estimated. 
That also supports the model of the n-distribution based on the 
return current. 

We expected that the bremsstrahlung from n-distribution 
would very likely manifest itself at low energies of the RHESSI 
spectral window. This was confirmed by the RHESSI fits of n- 
distribution discussed in the previous section. We note that in 
both cases the n-distribution significantly contributes to the X- 
ray spectrum in a rather narrow energy range, therefore the de- 
termined parameters of the n-distribution have large uncertain- 
ties. In all spectra the thermal component had to be considered 
to account for the emission in the ~ 6-10 keV energy range. 
Forcing the n-distribution component to higher energies instead 
of using the thermal component did not work: fits consisting of 
an n-distribution, two Gaussian lines (as approximations of the 
line complexes at ~ 6.7 and 8 keV), and a power-law high-energy 
tail did not lead to acceptable results. Finally, we also fitted the 
RHESSI spectra using a double power law and an isothermal 
component, i.e. without the n-distribution, to examine a possi- 
bility that the low-energy emission is due to a flat power-law 
component extending down to keV energies. Although an ac- 
ceptable fit of ~ 1 .3 was obtained, the resulting low-energy 
power-law part was very steep; i. e., <5l ~ H with Sl > S, thus 
resembling the n-distribution component. 
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Fig. 9. 2002 October 4 flare. Notation is the same as in Fig. [4] 
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Neither RHESSI nor RESIK spectra are spatially resolved. 
The radiation produced by the thermal component and the n- 
distribution comes from the coronal plasma, while hard X-ray 
bremsstrahlung origins primarily at the footpoints (Fig.fT2l. The 
25-50 keV contour is the largest one with 3-6 keV and 6- 
12 keV countours in the middle (Fig. [12], right) or in its upper 
part (Fig. Q2] left). Even though the contours of the 3-6 keV 
and 6-12 keV channels overlay each other, we suppose that the 
n-distribution and the thermal one are localized in different sub- 
structures of flaring plasma. There are several reasons that led to 
this conclusion. 

First, T determined from the RHESSI thermal component is 
always higher than the values of t from the n-distribution. This 
difference between log T and log r is usually 0.1- 0.2 dex, and it 
cannot be explained by the measured uncertainties (Figs. |4l fTJ [9] 
bottom left). The enhanced intensities of Si xnd lines observed 
by RESIK are the second argument for a different spatial lo- 
cation of the Maxwell and the n-distribution. The high ratio of 
Si xnd to Si xm lines needs a strong gradient of particles in the 
energy range above approximately 1.5 keV (Figs.[T]|2). The high 
number of electrons with the energy corresponding to the exci- 
tation energy of doubly excited state results in the high inten- 
sities of the satellite lines. On the contrary, the low number of 
electrons with energies higher than 2 keV results in the lower 



excitation rates, hence the lower intensities of the allowed lines. 
As a result, the observed line ratio of high n can be obtained. 
However, the relation between the intensities of the satellite and 
the allowed lines is not so straightforward because electrons with 
the energy above 2 keV influence also the ionization equilibrium 
and control the Sixm/Sixiv ratio. The n-distribution alone ful- 
fils the condition of a strong gradient of the electron distribution 
function in its high-energy range. 

The electron distribution composed of the n-distribution and 
the Maxwell one does not have such a strong gradient in its high- 
energy part because it has a much higher number of electrons 
there. Such a composed distribution has the same number of 
electrons with the energy corresponding the excitation energy 
of doubly excited state and therefore the same intensities of the 
satellite lines. But, it also has much higher intensities of the al- 
lowed lines due to the high number of particles with the energy 
higher than 2 keV. This must result in a much higher Si xm to 
Si xnd ratio than is observed. The line emission from different 
locations in the plasma is proportional to their emission mea- 
sures. The RHESSI spectra indicated that the emission measure 
of the n-distribution is much higher than for the Maxwell one 
(Tables lA?7llA.8l ): so the line spectrum of the n-distribution dom- 
inates o ver the Maxwell spectrum , which agrees with the obser- 
vations (Dzifcakov a et alJl201 lbl) . Finally, the electron densities 
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Fig. 11. RHESSI fits (from 4 keV) of the single power law (thick-target model), the thermal component, and the n-distribution 
(thin-target model) in the low energy range, together with residuals: the 2003 February 22 flare (left) and the 2003 January 7 flare 
(right). 



4201 



410| 

IT 
o 
in 
o 

i5 400 1 

>- 
o 



5 3901 
o 



380 



Reconstructed RHESSI Image 



22-Feb-2003 09:28:12.000 to 09:28:20.000 
Detectors: 3F 4F 5F 6F 7F 8F 
Energy Range: 12.0 - 25.0 keV 
Clean Total counts: 2.53E+04 




370 1 

50 



r- .o 

><D . 



60 70 80 90 
Heliocentric X (arcsec) 



CD 



IT) (£><£> 

in co 
cm 1 

i CD 

C CD 
CD—C 

-o-a 

«»» 

* CO ■■ 
■TSLU 



100 



o 

CD 
CO 

o 



-1601 



-170 



-180 



>- 

.2 -190 1 

C 

CD 
O 

g -200 1 

CD 
X 



Reconstructed RHESSI Image 



7-Jan-2003 23:28332.000 to 23:28:52.000 
Detectors: 3F 4F 5F 6F 7F 8F 
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Fig. 12. RHESSI images in the 12-25 keV channel with contours of 70% from channels: 3-6 keV (full line), 6-12 keV (dash-dot 
line), and 25-50 keV (long dash line): the 2003 February 22 flare (left) and the 2003 January 7 flare (right). The thin solid line 
marks the solar limb for the 2003 January 7 flare. 
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of the power-law distribution, as well as of the electrons pro- 
ducing the type III bursts, are too low compared to the other 
two assumed distributions to have any significant effects on the 
studied RESIK Si line spectrum; see also the recent results on 
a combined dis tribution composed of th e n-distribution and the 
power-law tail (Dzifcakov a et al.l l 201 lab . 

In summary, those facts point to different spatial origins for 
the n-distribution and the thermal emission. 

7. Conclusions 

In this analysis we have determined that during flares the plasma 
affected by an injected beam of accelerated electrons can have 
the non-thermal component in the 2-6 ke V region and that it 
can be described by the n-distribution. We derived the param- 
eters n and t of this n-distribution using the flare spectra from 
two satellite spectrometers: RESIK and RHESSI. The spectra 
were analysed by two distinct diagnostic methods. The soft X- 
ray lines of Si xiv, xm, xnd observed by RESIK provided us with 
the best possibility of analysing the behaviour of parameter n. 
The RHESSI analysis of the n-distribution is based on studying 
its continuous, bremsstrahlung radiation. Although the parame- 
ters n and r obtained from the RHESSI fits have large uncertain- 
ties, their values still overlap with the RESIK results quite well. 
The RESIK and RHESSI results indicate that departure from 
the Maxwell distribution at 2 - 6 keV energies can really occur. 
The non-thermal n-distribution and the power-law component 
are observed in good temporal correlation with the radio bursts 
detected in the 0.61 - 15.40 GHz range by RSTN. The thermal 
component is still present in RHESSI spectra. RHESSI images 
show that the emission from both thermal and non-thermal n- 
distributions comes from nearly the same area. Despite this we 
argue that in reality each distribution occupies different substruc- 
tures within the depicted area. Such a conclusion is supported 
by the character of the observed RESIK line spectra with high 
fluxes of satellite lines. Finally, our joint analysis of radio and 
X-ray emission also supports the model of formation of the n- 
distribution because of the return-current. 
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Appendix A: Complementary figures and tables 
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Table A.l. Time-weighted RESIK spectra for the C5.8 flare on February 22, 2003 and the determined ranges of the parameters n 
and log t. 



RESIK spectrum Time interval 
inUT 




n 






log T 
[K] 




GOES phase 





09:27:43- 


-09:28:31 


4.5 


-11 


-12 


7.13- 


-7.25- 


-7.26 


rise 


1 


09:28:31- 


-09:29:03 


3.1 


-5.2- 


-12 


7.07- 


-7.14- 


-7.27 


rise 


2 


09:29:03- 


-09:30:07 


2.5 


-3.0 


-12 


7.07- 


-7.10 


-7.14 


maximum 


3 


09:30:07- 


-09:31:11 


2.1- 


2.6 


-3.3 


7.02- 


-7.05- 


-7.08 


decay 


4 


09:31:11- 


-09:32:23 


1.9 


-2.5- 


-3.3 


6.94- 


6.99 


-7.03 


end 


5 


09:32:23- 


-09:33:51 


1.5- 


-2.5- 


-3.1 


6.89- 


6.94 


-6.99 


post-flare 


6 


09:33:51- 


-09:35:15 


1.5- 


-2.4 


-4.0 


6.87- 


6.94 


-7.03 


post-flare 


7 


09:35:15- 


-09:37:27 


2.2- 


3.4 


-9.5 


6.88- 


6.97 


-7.15 


post-flare 



Notes. Parameter values of n and r represent the determined values (bold) and their ranges. 



Table A.2. Results from RHESSI analysis of the 2003 February 22 flare. 



Time interval 


RESIK spectrum* 


Attenuator 




Thermal parameters 


Non-thermal parameters 


inUT 








EM 


kT 


io g r 


F \ 


6 E c 










[10 47 cm" 3 ] 


[keV] 


[K] 


[10 35 s-'] 


[keV] 


09:28:12-09:28:20 





AO 


0.97 


0.84 


1.88 


7.34 


0.29 


2.9 13 


09:28:28-09:28:36 


0, 1 


Al 


0.90 


2.9 


1.87 


7.34 


1.1 


5.4 21 


09:28:36-09:28:44 


1 


Al 


0.62 


6.0 


1.78 


7.31 


2.2 


5.7 19 


09:28:44-09:28:52 


1 


Al 


1.02 


10.5 


1.77 


7.31 


3.2 


6.9 19 



Notes. 



(*) 



See Table fAT 



Table A.3. Time-weighted RESIK spectra for the M4.9 flare on January 7, 2003. Notation is the same as in Table I A.l 



RESIK spectrum 


Time interval 
inUT 




n 






log T 
[K] 




GOES phase 





23:13:13-23:26:33 


1.0- 


3.1- 


8.1 


6.75 


-6.90- 


7.07 


pre-flare 


1 


23:26:33-23:28:25 


2.5- 


4.2- 


9.6 


6.87 


-6.98- 


7.12 


flare beginning 


2 


23:28:25-23:29:17 


3.3- 


6.5 


-12 


7.04 


-7.15- 


7.23 


rise 


3 


23:29:17-23:30:03 


2.9- 


3.8- 


6.2 


7.09 


-7.13- 


7.19 


rise 


4 


23:30:03-23:31:03 


4.6- 


6.1 


-12 


7.16 


-7.20- 


7.29 


rise 


5 


23:31:03-23:31:51 


7.7- 


11 


12 


7.25 


-7.29- 


7.30 


maximum 


6 


23:36:07-23:36:23 


11- 


11 


12 


7.22 


-7.22- 


-7.23 


maximum 


7 


23:40:39-23:41:03 


4.8- 


-7.3- 


-12 


7.16 


-7.22- 


7.28 


decay 


8 


23:41:03-23:42:03 


3.7- 


4.6 


6.0 


7.14 


-7.17- 


7.20 


decay 


9 


23:42:03-23:42:39 


2.9- 


3.5- 


4.9 


7.11 


-7.13- 


7.17 


decay 


10 


00:21:07-00:21:55 


1.6- 


2.0- 


2.6 


6.91 


-6.95- 


6.98 


decay 


11 


00:29:41-00:30:03 


1.2- 


1.9- 


2.8 


6.82 


-6.88- 


6.93 


decay 


12 


00:30:03-00:31:09 


1.9- 


2.4- 


2.9 


6.87 


-6.91- 


6.94 


decay 


13 


00:31:09-00:32:23 


1.0- 


1.2- 


1.5 


6.80 


-6.81- 


6.89 


decay 


14 


00:32:23-00:33:15 


1.0- 


1.9 


2.5 


6.78 


-6.85- 


6.89 


decay 


15 


00:33:15-00:34:09 


1.0- 


2.8- 


-11 


6.75 


-6.93- 


7.11 


decay 
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Table A.4. Results from RHESSI analysis of the 2003 January 7 flare. Notation is the same as in Table I A. 2l 



Time interval RESIK Spectrum* Attenuator xl Thermal parameters Non-thermal parameters 



in UT 








EM 

r i n47 ^,™-3i 

[iu cm j 


kT 

[Ke v j 


lop T 


LIU s J 


§ 


[Ke v j 


23:28:32-23:28:52 


2 


AO 


1.07 


0.64 


1.96 


7.36 


3.4 


7.0 


16 


23:29:16-23:29:36 


3 


Al 


1.05 


4.7 


1.73 


7.30 


2.9 


7.0 


19 


23:29:36-23:30:00 


3 


Al 


0.58 


5.6 


2.06 


7.38 


2.7 


6.9 


24 


23:30:00-23:30:28 


4 


Al 


0.89 


19. 


1.81 


7.32 


10.5 


6.3 


19 


23:30:28-23:31:00 


4 


Al 


0.94 


50. 


1.69 


7.29 


26. 


7.3 


18 


23:31:00-23:31:44 


5 


Al 


1.01 


67. 


1.89 


7.34 


27. 


7.1 


20 


23:32:00-23:32:32 




A3 


0.57 


45. 


2.38 


7.44 


2.9 


6.2 


26 


23:32:32-23:33:00 




A3 


0.84 


19. 


1.90 


7.34 


18. 


8.6 


22 



Notes. w See Table lA~3l 

Table A.5. Time-weighted RESIK spectra for the M4.0 flare on October, 4, 2002. Notation is the same as in Table IA.il 





Tirnp interval 

inUT 










log T 

[K] 




GDFS nhase 





05:35:19-05:36:03 


3.5 


-5.5- 


-12 


7.07 


-7.14- 


7.26 


rise 


i 


05:36:03-05:36:35 


2.1- 


-2.7- 


3.8 


6.97 


-7.01- 


7.06 


rise 


2 


05:36:35-05:37:01 


11 


11 


12 


7.21 


7.21- 


7.22 


rise 


3 


05:37:01-05:37:31 


3.1 


7.6 


-12 


7.09 


-7.22- 


7.28 


rise 


4 


05:37:31-05:38:01 


2.4 


-4.4- 


-12 


7.05 


-7.14- 


7.29 


maximum 


5 


05:38:01-05:38:31 


11 


11 


12 


7.24 


-7.24- 


7.25 


maximum 


6 


05:38:31-05:39:01 


11 


11 


12 


7.23 


-7.23- 


7.24 


decay 


7 


05:39:01-05:39:31 


2.6 


-4.2- 


-12 


7.09 


-7.16- 


7.29 


decay 


8 


05:39:31-05:40:01 


4.2 


-9.3- 


-12 


7.12 


-7.23- 


7.28 


decay 


9 


05:40:01-05:41:00 


2.5- 


-3.1- 


4.0 


7.07 


-7.10- 


7.14 


decay 


10 


05:41:00-05:42:00 


4.7 


-6.5 


-12 


7.16 


-7.20- 


7.29 


end 


11 


05:42:00-05:43:00 


3.3 


-4.1- 


5.9 


7.09 


-7.12- 


7.17 


post-flare 


12 


05:43:00 - 05:44:00 


2.8- 


-3.3- 


4.1 


7.08 


-7.10- 


7.13 


post-flare 


13 


05:44:00-05:45:00 


1.9 


2.1- 


2.8 


7.06 


-7.08- 


7.11 


post-flare 


14 


05:45:00 - 05:46:00 


2.7- 


3.1- 


3.8 


7.06 


-7.08- 


7.11 


post-flare 


15 


05:46:00-05:47:02 


2.6- 


3.0- 


3.6 


7.03 


-7.05- 


7.08 


post-flare 


16 


05:47:02-05:48:02 


1.8- 


-2.2- 


2.6 


6.98 


-7.00- 


7.02 


post-flare 


17 


05:48:02-05:49:02 


2.6- 


-2.9- 


3.5 


7.00 


-7.03- 


7.05 


post-flare 


18 


05:49:02 - 05:50:02 


3.0- 


-3.5- 


5.0 


7.03 


-7.06- 


7.11 


post-flare 


19 


05:50:02-05:51:06 


1.2 


-1.6- 


2.0 


6.93 


-6.96- 


6.98 


post-flare 


20 


05:51:06-05:52:02 


1.7 


-2.1- 


2.5 


6.92 


-6.95- 


6.98 


post-flare 


21 


05:52:02-05:53:06 


1.9- 


-2.3- 


-2.7 


6.92 


-6.95- 


6.98 


post-flare 


22 


05:53:06-05:54:06 


2.6- 


3.1- 


4.4 


6.99 


-7.02- 


7.08 


post-flare 
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Table A.6. Results from RHESSI analysis of the 2002 October 4 flare. Notation is the same as in Table lA.2l 



Time interval 


RESIK spectrum* 


Attenuator 




Thermal parameters 


Non-thermal parameters 


inUT 








EM 


kT 


logT 


F T 


S E c 










[10 47 cm" 3 ] 


[keV] 


[K] 


[10 35 s" 1 ] 


[keV] 


05:35:28-05:36:00 





Al 


1.18 


9.6 


1.84 


7.33 


586. 


5.8 8 


05:36:00-05:36:16 


1 


Al 


1.13 


51. 


1.89 


7.34 


27.1 


7.6 19 


05:36:32-05:37:00 


2 


A3 


0.88 


63. 


1.92 


7.35 


10.3 


7.9 22 


05:37:00-05:37:32 


3 


A3 


0.87 


105. 


1.86 


7.33 


10.5 


10.4 23 



Notes. w See Table lA~5l 



Table A.7. Details of the RHESSI spectrum fitting (2003 February 22, 09:28:12-09:28:20 UT) and comparison with the corre- 
sponding RESIK results (09:27:43-09:28:31 UT). 



Model 


Parameter 




Parameter values 




Energy range [keV] 


6-100 


4-100 


4-100 






0.97 


7.33 


1.12 


vth 


EM [10 47 cm" 3 ] 


0.84 


2.2 


0.70-0.85-0.90 


vth 


kT [keV] 


1.88 


1.45 


1.83-1.86-1.91 


photon_thick 


F T [10 35 electrons.s-'] 


0.29 


1.2 


0.09" -0.47 


photon_thick 


6 


2.9 


3.2 


2.8-2.9-3.0 


photon_thick 


E c [keV] 


13 


8.9 


9.8-32* 


thin_ndist 


EM„ [10 47 cm" 3 ] 






5.0" -5.1 


thin_ndist 


n 






2-9 (ll c ) 


thin_ndist 


kr [keV] 






1.0-1.5 (1.63 r ) 


RESIK 


n 






4.5-11-12 


RESIK 


kr [keV] 






1.16-1.53-1.57 



Notes. Columns three to five list parameter values of three types of RHESSI fits (original, i.e. vth + photon-thick, original within the enlarged 
energy range, and the one combined with the n-distribution), see Sect. [5] and left panels of Figs. IA. ill A.2I andfTT] Last column shows the fitted 
parameters (bold), together with their limits. <a) lower limit only (b) upper limit only w values outside estimated lcr range 



Table A.8. Details of the RHESSI spectrum fitting (2003 January 7, 23:28:32-23:28:52 UT) and comparison with the correspond- 
ing RESIK results (23:28:25-23:29:17 UT). Notation is the same as in Table |A71 



Model 


Parameter 




Parameter values 




Energy range [keV] 


6-40 


4-40 


4-40 




xl 


1.07 


5.38 


1.31 


vth 


EM [10 47 cm- 3 ] 


0.64 


4.1 


0.3-0.35-0.41 


vth 


kT [keV] 


1.96 


1.18 


2.3-2.44-2.7 


photon.thick 


F T [10 35 electrons.s-'] 


3.4 


9.6 


1.6-2.4-3.1 


photon.thick 


<5 


7.0 


7.2 


5.8-6.5-6.7 


photon.thick 


Fc [keV] 


16 


14 


15-15.5-16 


thin_ndist 


EM B [10 47 cm" 3 ] 






0.5° -4.5 


thin_ndist 


n 






5.0° -6.5 


thin_ndist 


kr [keV] 






1.5-1.72-2.6 


RESIK 


n 






3.3-6.5-12 


RESIK 


kT [keV] 






0.94-1.22-1.46 



Notes. See also Sect.[5]and right panels of Figs. lA. HlA.21 andfTT] 
(fl) lower limit only 




Fig. A.l. RHESSI fits (from 6 keV) of the single power law (thick-target model) and the thermal component together with residuals: 
the 2003 February 22 flare (left) and the 2003 January 7 flare (right). 
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SPEX HESSI Count Flux vs Energy with Fit Function, Interval SPEX HESSI Count Flux vs Energy with Fit Function, Interval 



i Detectors: 4F 

22-Feb-2003 09:28:12.000 to 09:28:20.000 (Data with Bk) — 
22-Feb-2003 09:28:12.000 to 09:28:20.000 (Bk) — 




Fit Interval !0 Chi-sq u are = 7.33 

— vth+phbton_thicK+pileup_mod+drm_mod 
vth 0.0222,1.45,1.00 full chian 1 .26e-04 
photon thick 1.16,3.20.200. .18.0,8.88. 1.00e+( 

— pileupimod 1.19,0.900,1.00,10.0,3.00,0.500 
drm_rrtod 0.608,0.0777,1.00 



i Detectors: 4F " 

7-Jan-2003 23:28:32i000 to 23:28:52.000 (Data with Bk) — 
7-Jan-2003 23:28:32.000 to 23:28:52.000 (Bk) — 




Fit Interval !0 Chi- square = 5.38 

— vth+phbton_thick+drm_mod 
vth 0.0403,1.18,1.00 full chian 1 .26e-04 
photorLthick 9.62,7.24,200.,8.00,1 4.0,1 .00e+03 

— drm_rrtod 0.837,0.00851,1.00 



Energy (keV) 
Residuals for Interval 



Energy (keV) 
Residuals for Interval 




Fig. A.2. Widened RHESSI fits (from 4 keV). Notation is the same is in Fig. IA.ll 



